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1. Introduction. 

The separation of variables is one of the most universal methods for solving com- 
pletely integrable (classical and quantum) models. It has been applied successfully 
to the study of a large number of finite- dimensional integrable Hamiltonian systems 
(FDIHSs) (see, for example, [1-12]), as well as infinite dimensional integrable Hamil- 
tonian systems in the determination of finite-dimensional quasi-periodic solutions (see, 
for example, [13-18]). In many cases the separation of variables of integrable classical 
systems prepares the passage to the corresponding quantum systems. For the classical 
integrable systems subject to the inverse scattering method, the standard construction 
of the action-angle variables using the poles of the Baker-Akhiezer function is in fact 
equivalent to the separation of variables [4] . 

For a FDIHS, let m denote the number of degrees of freedom, and P^, i = 1, m, be 
functionally independent integrals of motion in involution, the separation of variables 
means to construct m pairs of canonical separated variables v^, Ufe, /c = 1, m, [2,3,4] 

{uk.ui} = {vk,vi} = {vk.ui} = Sku k,l = l,...,m, (1.1) 

and m functions fk such that 

fk{uk,Vk,Pi,...,Pm) = 0, k=l,...,m. (1.2) 

The equations (1.2) are called separated equations, which give rise to an explicit fac- 
torization of the Liouville tori. 

For the FDIHSs with the Lax matrices admitting the r-matrices of the XXX, XX Z 
and XYZ type, there is a general approach to introduce canonical separated variables 
[2,3,4,8]. The corresponding separated equations enable us to express the generating 
function of canonical transformation in completely separated form as an abelian integral 
on the associated invariant spectral curve. The resulted linearizing map is essentially 
the Abel map to the Jacobi variety of the spectral curve, thus providing a link, through 
purely Hamiltonian methods, with the algebro-geometric linearization methods given 
by [19-22]. 

An important feature of the separation of variables for a FDIHS is that the number of 
canonical separated variables Uk should be equal to the number m of degrees of freedom. 
In some cases, the number of resulted by the normal method may be less than m 
and so some additional canonical separated variables should be introduced. So far very 
few models in these cases have been studied. These cases remain to be a challenging 
problem [4] . In recent years binary constrained flows of soliton hierarchies have attracted 
attention (see, for example, [23-29]), whose basic idea was described in [30]. The binary 
constrained flows are a kind of FDIHSs for which the method presented in [2,3,4,8] is 
not valid. The degree of freedom for binary constrained flows admitting 2x2 Lax 
matrices is an even natural number usually denoted by 2N . The method in [2,3,4,8] 
allows us to introduce only N pairs of canonical separated variables ui,...,un and 
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vi, ...,vn via the Lax matrices. In this paper we propose a new method for determining 
additional N pairs of canonical separated variables and separated equations for binary- 
constrained flows. The main idea is to construct two functions B{X) and A{X) defining 

un+1, ■■■,U2n by the set of zeros of B{\) and v^^^k = ^iuN+k)- To keep the canonical 
conditions (1.1) and the requirement for the separated equations (1.2), it is found that 

certain commutator relations should be imposed on -S(A), A(A) and A[\) has some link 
with the generating function of integrals of motion of binary constrained flows, which 
provides a way to construct the B{\) and ^(A). In fact, we have to modify the original 
approach for introducing lii, ...,U7v and vi,...,vn so that ui, ...,U2n and vi, ...,V2n are 
canonical conjugated. Having produced the separation of variables, we further construct 
the Jacobi inversion problems for binary constrained flows. This method is somewhat 
different from that for introducing canonical variables presented in [31] and can be 
applied to more binary constrained flows. 

Briefly, separation of variables can be characterized as a reduction of a multidimen- 
sional problem to a set of one-dimensional ones. The separation of variables of soliton 
equations in this paper contains two steps of separation of variables. The first step is to 
factorize (1 + 1)— dimensional soliton equations into two commuting x— and t— FDIHSs 
via binary constrained flows, namely the x— and t— dependences of the soliton equations 
are separated by the x— and t— FDIHSs obtained from the x- and t-binary constrained 
flows. The second step is to produce separation of variables for the x— and t— FDIHDs 
by our method to be proposed later on. Finally, combining the factorization of soliton 
equations with the Jacobi inversion problems for x— and t— FDIHSs enables us to estab- 
lish the Jacobi inversion problems for soliton equations. We will present the separation 
of variables for the KdV hierarchy, the AKNS hierarchy and the Kaup- Newell hierarchy 
via their binary constrained flows. In fact, we employ our method in a little different 
way for those three cases. 

In section 2, we recall the binary constrained flows of KdV hierarchy and present fac- 
torization of the KdV equations into the x- and t-binary constrained flows. By means 

of the Lax matrix M = ( for the binary constrained flows, the method 

\C{X) -A{X)J 

in [2,3,4,8] allows us to define only pairs of canonical variables ui,...,un by the set 
of zeros of -B(A) and Vk = 2A{uk). We propose a new method to construct two new 
functions -B(A) and A{X) for introducing un+1i ■■■■,U2n by the set of zeros of B{X) and 
VN+k = A{uN+k)- The construction of B{X) and A{X) is based on an observation that 
the canonical conditions (1.1) need certain commutator relations between A(A),S(A), 
and the requirement for the separated equations (1.2) links A{X) with another gener- 
ating function of integrals of motion. To guarantee that vi, ...,V2n and ui, ...,U2n are 
canonical conjugated, we also have to modify the original way for introducing ui, .., un 
and vi, ...,vn- Then we establish the Jacobi inversion problems for the x- and t-binary 
constrained flows. Finally, these Jacobi inversion problems together with the factor- 
ization of the KdV equations give rise to the Jacobi inversion problems for the KdV 
equations. In section 3, the factorization of the AKNS equations is given. Since B{X) 
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for the binary constrained AKNS flows, unlike the -B(A) for the binary constrained KdV 
flows, has only N — 1 zeros, we have to modify the method proposed in section 2 in 
order to find 2N pairs of canonical variables for the binary constrained AKNS flows. In 
section 4, we present the factorization of the Kaup-Newell equations. Since the com- 
mutator relations of ^(A), -B(A) and C(A) for the binary constrained Kaup-Newell flows 
are quite different from those for both the binary constrained KdV flows and the binary 
constrained AKNS flows, we need to further modify the method in sections 2 and 3 in 
order to flnd the separation of variables for the Kaup-Newell equations. Finally some 
remarks are made in section 5. 

2. Separation of variables for the KdV equations. 

In this section, we use the binary constrained flows of KdV hierarchy to illustrate our 
method of introducing canonical separated variables. Then we show how to produce 
the separation of variables for the KdV equations. To make the paper self-contained, 
we flrst briefly describe the binary constrained flows of the KdV hierarchy [26] . 

2.1 Binary constrained flows of the KdV hierarchy. 

Let us start from the Schrodinger equation [32] 

(f>xx + (A -I- u)(f) = 
which can be rewritten as following spectral problem 

<P. = U{u,X)<P, C/(»,A)=(_/_^ J), ^=(^^). (2.1) 

Its adjoint representation reads 

= [U, V] = UV- VU. (2.2) 

Set 



oo , 

' X-\ (2.3) 



i=0 



Ci di 



Equation (2.2) yields 



aQ = bQ = 0, Co = -1, ai = 0, bi = 1, ci = —u, 
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and in general 



1 I. 1 1/ 2^ 

a2=^Ux, b2 = --u, C2 = -{Uxx + u ),..., 



bk+i = Lbk = -^^'^ > = -^bk,x, (2.4a) 



where 



Set 
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ck = ~h,xx - h+1 - hu, /c = l,2,---, (2.4b) 



i=0 

and take the time evolution law of as 



y(")(ii,A)0. (2.6) 



Then the compatibility condition of the equations (2.1) and (2.6) gives rise to the n-th 
KdV equation which can be written as the infinite-dimensional Hamiltonian system 

= -2bn+2,x = dL^u = d^, (2.7) 
where the Hamiltonian is given by 

The matrix V determined by (2.2) and (2.3) also satisfies the adjoint representation 
of (2.6) 

Vt^ = [V^^\V], (2.8) 

when u satisfies (2.7). 
For n = 1 we have 



0,, = fW(«,a)0, = ( 



(2.9) 



and the equation (2.7) for n = 1 is the well-known KdV equation 



Ut^ = -^{Uxxx + Quux). (2.10) 



The adjoint spectral problem reads 



V'. = -f/^(^,A)V, ^=(^M- (2-11) 



We have [26] 
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where (3 is some constant. 

The binary x-constrained flows of the KdV hierarchy (2.7) consist of the equations 
obtained from the spectral problem (2.1) and the adjoint spectral problem (2.11) for 

distinct real numbers A-,- and the restriction of the variational derivatives for the 
conserved quantities Hk^ (for any fixed ko) and A^-: 

^i,x = ^2, ^2,x = -A^i - u^i, (2.13a) 

^l,x^ M!2+U^2, ^2,x^-^l, (2.13b) 

^ - E ^ = -2^^0+2+ < >= 0. (2.13c) 

Such a constraint (2.13c) has been recognized as a symmetry constraint [25,26,30]. 
Hereafter we denote the inner product in by < ., . > and 

= (4i, • • • ,0iiv)'^, *i = (^ii, • • • , V'iiv)'^, ^ = 1,2, A = (im5f(Ai,- • • , Aat). 
For /co = 0, we have 

1 1 

h2 = --u= -<m2,^i>, i.e., 'u = -<*2,*i>. (2.14) 

By substituting (2.14) into (2.13a) and (2.13b), the first binary x-constrained fiow be- 
comes a finite-dimensional Hamiltonian system (FDHS) [26] 

dFi ^ dFi ^ dFi ^ dFi 

with the Hamiltonian 

Fi ^'i, $2 > - < $1 > +^ < *2, $1 >^ • 

The binary t^-constrained fiows of the KdV hierarchy (2.7) are defined by the replicas 
of (2.6) and its adjoint system for N distinct real number Aj 

(2.16) 

as well as the n-th KdV equation itself (2.7) in the case of the higher-order constraint for 
/cq > 1. Under the constraint (2.14) and the x-FDHS (2.15), the binary ti-constrained 
flow obtained from (2.16) with V'^'^^ given by (2.9) can also be written as a ti-FDHS 

dF-) dF-y dF-) dF-) 

^Mi Q^^^ ^2,ti ^i,t^ Q^^^ V^-^n 
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with the Hamiltonian 

F2 = - < A^*2, $1 > + < A*i, ^2 > +^ < *2, *i >< A*2, ^1 > 

1 1 

2 o 

The Lax representation for the x-FDHS (2.15) and the ti-FDHS (2.17) can be de- 
duced from the adjoint representation (2.2) and (2.8) by using the method in [33,34] 

M, = [^,M], Mt„ = [y("),M], (2.18) 

where JJ and V'^"') are obtained from JJ and V^'^^ by a substitution of (2.14), and the 
Lax matrix M is given by 

j = l J 

The equation (2.18) imphes that \TrM'^{X) = ^^(A) + B{X)C{X) is the generating 
function of integrals of motion for (2.15) and (2.17). A straightforward calculation 
yields 

N p p2 

A\X) + B{X)C{X) ^ P(A) = -A + E[ A^A- + (A^l^]' ^^.20) 

where Pj,j = 1, ...,2iV, are 2N independent integrals of motion for the FDHSs (2.15) 
and (2.17) 

1 x 1 

+ 8 X] \ . _ \ [(^IJ^lj ~ V'2i02j)(V'lfe01fc - -02/002^) +4V'lj02jV'2/c01fc], 



j = l,...,N 

(2.21a) 

Pn+j = ^{^lAij + V^2i02i), j = 1, N. (2.21b) 



It is easy to verify that 

N N 



Fi = 2^P,-, F2 = 2^(A,P, + P^+,.)- (2.22) 

3=1 3=1 



8 



YUNBO ZENG AND WEN-XIU MA 



With respect to the standard Poisson bracket it is found that 

{A(A), = {S(A), Bifi)} = {C(A), Ciii)} = 0, (2.23a) 

{A(A),S(/.)}= 2p^[S(/.)-S(A)], (2.23b) 

{A{X),C{f,)} = ^^^^[C{X) - C{fi)], (2.23c) 

{BiX),Ci^i)}=-^[Ai^^)-AiX)]. (2.23d) 
A — ^ 

It follows from (2.23) that 

{A\X) + S(A)C(A), A^fi) + B{fi)C{fi)} = 0, 

which implies that Pj,j = 1, ...,2N, are in involution: 

{Pk,Pi} = 0, kJ = l,...,2N. 

Therefore the FDHSs (2.15) and (2.17) are integrable and commute with each other. 
The construction of (2.15) and (2.17) ensures that if (\l/i, \l/2, ^i, ^2) satisfies the finite- 
dimensional integrable Hamiltonian systems (FDIHSs) (2.15) and (2.17) simultaneously, 
then u defined by (2.14) solves the KdV equation (2.10). 

In general, by substituting (2.14) and using (2.15), the t„-constrained flow (2.16) 
becomes a t„-FDIHS and the n-th KdV equation (2.7) is factorized by the x-FDIHS 
(2.15) and the ^FDIHS. Set 



A^{X) + B{X)C{X) = Xj2Fk>^~^, (2-24a) 

fc=0 

where Ffc, A; = 1, 2, are also integrals of motion for both the a;-FDHSs (2.15) and the 
t„-binary constrained flows (2.16). Comparing (2.24a) with (2.20), one gets 

N 

Fo = -l, Fi=0, Ffc = ^[A)-2p, + (A;-2)Aj-3p2^^.], A; = 2,3,.... (2.24b) 

By employing the method in [34,35], the i^-FDIHS obtained from the i^-constrained 
flow (2.16) is found to be of the form 

-,- dFn+i dFn+i dFn+i dFn+i . . 
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with the Hamiltonian 



m=0 li + ...+lm+i=n+2 

where li > 1, Im+i > 1, cto = 1? cti = |) ck2 = and [34,35] 

m — 2 ^ m — l 

am = '2am-i + ^ q;/q;,„_«_i ^ ^ m > 3. (2.25c) 



2 

/=i 1=1 

The n-th KdV equation (2.7) is factorized by the x-FDIHS (2.15) and the tn-FDIHS 
(2.25). 

For example, for the second equation in the KdV hierarchy (2.7) with n = 2 

Ut2 = Y^iuxxxx + iOuu^x + 5ul + 10u^)x, (2.26) 

the Hamiltonian F3 for the t2-FDIHS reads 

~ 1~ ^ 1 ^ 

F3 = 2F4 + -Fi = 25^(AfP, + 2A,P^^,.) + -Q^Pjf- (2-27) 

Then the second KdV equation (2.26) is factorized by the x-FDIHS (2.15) and the 
i2-FDIHS with the Hamiltonian F3. 

2.2 The separation of variables for the KdV equations. 

An effective way to introduce the separated variables Vk,Uk and to obtain the sep- 
arated equations is to use the Lax matrix M and the generating function of integrals 
of motion. For the FDIHSs (2.15) and (2.17), we can define the first N pairs of the 
canonical variables Uk,Vk,k = 1, ...,N, by the method [2,3,4,8]. The commutator rela- 
tions (2.23) and the generating function of integrals of motion (2.20) enable us to define 
ui,...,un by the set of zeros of B{X) 

where 

N N 

R{X) = 1[{X- Uk), K{X) = l[iX- Xk), 

k=l k=l 
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and vi, vn by 

Vk = Ai{uk), k = l,...,N, (2.28b) 

where 

MX) = 2A(A) = i V ^i.'^i.- ^2.'^2. . 

As we will see later, the commutator relations (2.23) guarantee that ui,...jUn and 
Vi,...,Vn satisfy the canonical conditions (1.1). Then substituting Uf- into (2.20) gives 
rise to the separated equations 



Vk 



= Ai(uk) = 2y/P(uk) = 2 



N 



X, 



+ 



p2 



k = 1,...,N. 



Now the reason for taking our choice of B{X) and Ai{X) becomes apparent. 

The FDIHSs (2.15) and (2.17) have 2N degrees of freedom, therefore we need to 
introduce the other N pairs of canonical variables Vk,Uk.,k = N + 1, 2N. The main 
idea is to construct two suitable functions B{X),A{X) in order to define un+i, ■■■,U2n 
by the set of zeros of B{X) and vn+i, V2n by v^+k = A{uN-irk)- The above way for 
introducing Uk,Vk, k = 1, N stimulates us to impose two requirements on B{X) and 
A{X) in order to construct them. First, the canonical conditions (1.1) require that B{X) 
and A{X) satisfy 



{5(A), B(^)} = {B{X),B{f^)} = {A{X),A{f^)} = {A{X),B{^)} = {A{X),AM} = 0, 

(2.29a) 

{I(A),B(/i)} = ^[S(/x) -5(A)], (2.29b) 

A — fJ, 

{Ai{X),B{fi)} = 0. (2.29c) 

The second requirement is that the equation VN+k — A{uN+k) should give rise to the 
separated equations. Notice that Pn+j given by (2.21b) are integrals of motion for 
the FDIHSs (2.15) and (2.17), we can construct another generating function A{X) of 
integrals of motion by 

X^^^ ^ 1 ^ ^i.0i. +^^2,02, ^ 2 ^2 (2.30a) 
3=1 ^ i=i ^ 



We may use A{X) to define v^j^i, ...,V2n since substituting u^+k into the equation 
(2.30a) immediately leads to the separated equations for VN+k and wjv+fc- It is easy to 
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see that {A{X),B{^)} = {A(A),A(/x)} = {A{X),Ai{^)} = 0. We look for B{X) in the 
form 



2 ^ A- A 

By requiring B{X) to satisfy (2.29a) and (2.29b), one gets Si — 1, S2 — S3 — 0, i.e. 



^W = i + ^Ei3t:- (2.30b) 



2 ^ A - A, 
But B{X) doesn't fit (2.29c). In fact, we have 

{A^{X),B{^i)} = -^[B{^i)-B{X)]. (2.31) 

A — jJL 

However, (2.29a), (2.29b) and (2.31) enable us to replace Ai(A) by A(A) 

A(A)^^,(A)-I(A) = -f;^, (2.32) 

j=i 

and we will redefine Vk by Vk = A{uk)- 

Then a straightforward calculation shows that -B(A) = B{X), A{X), B{X), A{X) satisfy 
the following required commutator relations 

{B{X),B{fi)} = {BiX),Bifi)} = {A{X),A{fx)} = {1(A), !(//)} = 0, (2.33a) 

{B{X),B{i^)} = {B{X),A{i^)} = {B{X),A{i^)} = {A(A),I(ai)} = 0, (2.33b) 

{AiX),Bifi)} = -^[B{f,)-BiX)], {I(A),5(/.)}=-^[S(/.)-S(A)]. (2.33c) 

A — fl A — /J, 

We have the following proposition. 

Proposition 1. Assume that Xj , (pij , ipij G R, i = l,2,j = 1,...,N. Introduce the 
separated variables ui, ...,U2n by the set of zeros of B{X) and B{X): 



i^-Bi^^.l^^^m, (,34a) 



5W = l + iEx^^|^, (2.34b) 
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with 

N N 

R{x) = - ^^)' ^(^) = - ^^+'^)' 

k=l k=l 

and vi, ...,V2N by 

Vk = A{uk) = Ai{uk) - A{uk) = - V in^, k^l, N, (2.34c) 
^^^^+.^A(n^W-|E ^^^^^^^^^^f^^ /c = l,...,iV. (2.34d) 

^ ~{ UN+k - Aj 

If ui, ...,un, are single zeros of B{X), then vi,...,V2n and ui,...,U2n are canonically 
conjugated, i.e., they satisfy (1.1). 



Proof. Notice that 



lim -B(A) = — oo, lim 5(A) = oo, 



it is easy to see that 

un+1 < Ai < un+2 < A2 < ... < U2N < Ajv- (2.35) 

We have B (uk) 7^ 0, B'{uj^^k) ^ 0. Hereafter the prime denotes the differentiation 
with respect to A. In what foUows, we take A;, / = 1, A^. It foUows from (2.33b) that 

= {Uk,B{uN+l)} = B'{uN+l){Uk,UN+l} + {Uk,B{n)}\fj,=uj^^^, 
= {B{uk),UN+l} = B'{uk){uk,UN+l} + {B{X),UN+l}\x=Uk, 

= {B{uk),B{uN+i)} = B'{uk)B'{uN+i){uk, un+i} + B'{uk){uk, B{ii)}\ij,=un+i 
+B'{uN+l){B{X),UM+l}\ + {B{X),B{fi)}\ 

= B (uk)B'(uN+i){uk, UN+i} + B {uk){uk, B{ii)}\ 

IJ,=Un+1 

+B'{un+i){B{X) , UN+l } I \=uk , 

which together lead to {uk,UN-^i} = 0. Similarly, {uk,ui} — 0, {ujv+fe? '^^JV+i} = 0. 
Using (2.33b), (2.33c) and above results, one gets 

{vk,B{ii)} \i^=ui= {A{uk),B{ii)} \f,^ui 

= A'{uk){uk,B{ii)} +[{A{X),B{fj,)} \x=uk] \i^=ui, 
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= A {uk)[{uk, B{ui)} - B {ui){uk,ui}] + [{A{\),B{^j)} 



-Uk\ \tJ' = Ul 



[{A{X),B{t,)} U.„J = ^^^^ ^^"""^ \,=u,= -SkiB'M, 

Uk - II 



and 
then 



0^{vk, B{ui)} = B {ui){vk,ui} + {vk, B{n)} 



{vk.ui} = {vk,B{ij)) \^=ui= hi^Jj^ = Ski- 

B [ui) B (ui) 

In the same way, one gets {vjv+fc) i^JV+z} = ^ki- The foUowing equahties 

{vk,UN+l} ^ {A{Uk),UN+l} = {A{X),Un+i} |A=Mte, 

= {A{X),B{uN+i)} = B'{uN+i){A{X),UN+i}, 

yield {vk, un+i} = and similarly {vN+k, n} = 0. 
Finally, 

{vk,VN+l} = {A{Uk),A{uN+l)} 



= A (uk){uk,A{ij)}\ 

H=UN+i +^ {uN-\-i){.-^{^)i'^N+i}\\=uu 
= A'{uk)[{Uk, VN+l} - A'{uN+l){Uk, UN+l}] 
+A'{uN+l)[{vk,UN+l} - A'{uk){uk,UN+l}] = 0, 

similarly 

{vk, vi} = A'{uk){uk, vi} + A'{ui){vk, ui} = -A'{uk)Ski + A'{ui)Ski 

and {vN+k,VN+i} — 0. This completes the proof. 
It follows from (2.34a) and (2.34b) that 



or 



^ ^^'^ ^i?(A,)K'(A,) 



Also (2.34a) results 

N 

= - < *2,^i >= 2j2{uj-Xj). 
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We now present the separated equations. By substituting Uk into (2.20), u^+k into 
(2.30a) and using (2.34), one gets the separated equations 



Vk = Ai{uk) - A{uk) = 2^P{uk) - A{uk) 



= 2. 



N p 



VN+k = A{uN+k) = 2 V ^ = 1' ^- (2-38b) 

^ UN+k - Xj 

Replacing Vk by the partial derivative of the generating function S of the canonical 
transformation and interpreting the Pi as integration constants, the equations (2.38) 
give rise to the Hamilton- Jacobi equations which are completely separable and may be 
integrated to give the completely separated solution 

S{uu U2n) = / (2/P(A) - A{\))d\ + / A{\)d\] 

k=l 



,/p(X)dX -TPm+M I ~ ^\ |]. (2.39) 

jH' UN+k - Aj 



fe=l " i=l 

The linearizing coordinates are then 



ds A n 1 



By using (2.22), the linear flow induced by (2.15) is then given by 



(2.40b) 



dFi dFi 
Qi = li + x—- ^ ^i + 2x, QN+i = '^lN+i + x— = 27iv+i, z = 1, AT. (2.41) 

Hereafter 'ji,i= 1, 2N, are arbitrary constants. Combining the equation (2.40) with 
the equation (2.41) leads to the Jacobi inversion problem for the FDIHS (2.15) 

E/ 7^ rw^^^ = ^^ + 2a;, z = l,...,iV, (2.42a) 

fc=W [X- Ai)^yP{A) 
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y ^M+i dX-ln\ "\ |]=7iv+», z=l,...,iV. (2.42b) 

The 'i/'2i and < \l/2,$i > defined by (2.36) and (2.37) are the symmetric functions 
of Uk,k = 1,...,2N. If, by using the Jacobi inversion technique [19], 0ij,'i/'2j and 
< ^2,^1 > can be obtained from (2.42), then ^2j, V'lj can be found from the first and the 
last equation in (2.15) by an algebraic calculation, respectively. The (<?!>ij, 02^? V'lj) V'2j) 
provides the solution to the FDIHS (2.15). 

By using (2.22), the linear flow induced by (2.17) is then given by 

Qi = li + ^^^1 = 7i + 2Aiti, 
dF 

Qn+z = 2jN+i + ^ ^1 = 27iv+i + 4Pjv+iti, i = l,...,N, (2.43) 

where 7i are arbitrary constants. Combining the equation (2.40) with the equation 
(2.43) leads to the Jacobi inversion problem for the FDIHS (2.17) 



k=l 

N 



El/ 

fc=i 



— ——==dX - In = 7iv+j + 2PN+iti, I = 1, N. 

(A - Xi)^^yP{X) UN+k - Aj 



(2.44b) 

Finally, since the KdV equation (2.10) is factorized by the FDIHSs (2.15) and (2.17), 
combining the equation (2.42) with the equation (2.44) and using (2.37) give rise to the 
Jacobi inversion problem for the KdV equation (2.10) 



J2 I _ ^ )/P(A) '^'^ - 7i + 2x + 2Ai^i, i = l,...,N, (2.45a) 



k=l 

N 



"^^^ dX-ln\ ^ '\ |] =7jv+^ + 2PJv+^tl, i = l,...,iV. 



{X-XiY./P{X) 'uN+k-Xi 

(2.45b) 

Notice that u defined by (2.37) is the symmetric function of tt/t, = 1, A^. If, by using 
the Jacobi inversion technique [19], u can be found in terms of Riemann theta functions 
by solving (2.45), then u provides the solution of the KdV equation (2.10). 

In general, since the n-th KdV equation (2.7) is factorized by the x-FDIHS (2.15) and 
the tr^-FDIHS (2.25), the above procedure can be applied to find the Jacobi inversion 
problem for the n-th KdV equation (2.7). We have the following proposition. 
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Proposition 2. The Jacobi inversion problem for the n-th KdV equation (2.7) is given 
by 



/ ; dX = 7i + 2a; 

(A-A,)a/P(A) 



+*nE(or"'«- E \'r'-^Fi,...Fi^, z=l,...,7V, (2.46a) 



•2 

m=0 ii + ...+«^+i=n+2 

<-Uu T) \ 



k=l 

n 



' (A - Ai)2i/P(A) UN+k-\ 



1 ^ ^ 



'2' 

m=0 /i+...+/^+i=n+2 



i=l,...,A^, (2.46b) 

where li > 1, > 1 and Fi^, ■■■Fi^, are given by (2.24b). 

For example, by using (2.27), the Jacobi inversion problem for the second KdV 
equation (2.26) is given by 

J2 dX = ^, + 2x+i2Xt + J2Pj)t2, i = l,...,N, (2.47a) 

t[J (A-A,)2/P(A) UM+k-X^^ 

(2.47b) 

The u solved from the Jacobi inversion problem (2.47) provides the solution for the 
second KdV equation (2.26). 

The Jacobi inversion problem for the KdV hierarchy in our case is somewhat different 
from that derived by means of the stationary equations of the KdV hierarchy [36] , since 
there is an additional term —ln\uk — Xi\ + ln\uN+k — Ai| in (2.46b). 

3. The separation of variables for the AKNS equations. 
3.1 Binary constrained flows of the AKNS hierarchy. 

For the AKNS spectral problem [37] 
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its adjoint representation (2.2) and (2.3) yield 

ao = -1, &o = Co = 0, ai = 0, h = q, ci = r, a2 = ^qr, 
and in general 

ak = d~^{qck - rbk), /c=l,2, 

T _ ^ ( d — 2rd~^q 2rd~^r \ 
2\ -2qd-^q -d + 2qd-^r J ' 



(3.2) 



Take 



0,„=y(-)(«,A)0. !/(-)(«, A) = f;f^^ })x-\ (3.3) 

Then the compatibility condition of equations (3.1) and (3.3) gives rise to the AKNS 
hierarchy 



where the Hamiltonian Hn and the Hamiltonian operator J are given by 

-0 0^)' fe)^^' 

For 77. = 2 we have 

and the AKNS equation (3.4) for n = 2 reads 

(lt2 = ~Qxx + q^r, = ^ra;x - r^q- (3.6) 

The adjoint AKNS spectral problem is the equation (2.11) with U given by (3.1). 
We have [25] 

^ ^ or^r/^'/'i^i <Pi^2 \ dU{u,X) ^ nf^i<p2\ . . 

The binary ar-constrained flows of the AKNS hierarchy (3.4) are defined by [25,29] 
$1,^ = -A$i + g$2, ^2,x = r$i + A$2, (3.8a) 
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= A*i -r*2, ^'2,x = - A*2, (3.8b) 

r'i:^^f?!'0-f;!'1?:^o. (3.80 



5u j^i ^'^ V^fco/ V<*2,$1> 



For fco = 1, we have 



^1/ V?/ V<*2,*i> 



= 0. (3.9) 



By substituting (3.9) into (3.8a) and (3.8b), the first binary a;-constrained flow becomes 
a x-FDHS [25] 

dFi ^ dFi ^ dFi ^ dFi 

with the Hamiltonian 

Fi =< A^'2, $2 > - < A*i, $1 > + < *2, $1 >< $2 > • 

Under the constraint (3.9) and the FDHS (3.10), the binary t2-constrained flow ob- 
tained from (3.3) with V^"^^ given by (3.5) and its adjoint equation for N distinct real 
number can also be written as a i2-FDHS 

dFo dFo dFo dFo 

with the Hamiltonian 

F2 =< A^^'2, $2 > - < A^^'i, $1 > + < ^'2, $1 >< A^'i, $2 > 

+ < A*2,$l >< *1,$2 > -^(< *2,$2 > - < >) < ^'2,$l >< ^'l,$2 > • 

The Lax representation for the FDHSs (3.10) and (3.11) which can also be deduced 
from the adjoint representation (2.2) and (2.8) are presented by (2.18) with the entries 
of the Lax matrix M given by [29] 

N 

A{X) = -1 + ^ E (3.12a) 

A A An 

3 = 1 ^ 



.7=1 .7=1 



(3.12b) 
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A straightforward calculation yields 



^ Pi 



A\X) + B{\)C{X) ^ P(A) = 1 + Etx^ + (A^S^^' ^^-^^^ 

where Pj,j = 1, ...,2N, are 2A/' independent integrals of motion for the FDHSs (3.10) 
and (3.11) 

+ j \ . 1 \ [fe01j-V'2j02j)(V^lfc^lfc-'02fe02fe)+4'01j02jV'2fe01fc], j = 1, 



2 f-f. Xj — Xk 



(3.14a) 

Pn+j = li-ipiAij + ip2j4>2j), j = 1, N. (3.14b) 



It is easy to verify that 

N , N 



Fi = 5^(A,P, + P^+^) - liY^Pjf. (3.15a) 



AT , AT AT -.A/" 



j=i j=i j=i j=i 

With respect to the standard Poisson bracket it is found that 

{A(A), A(/x)} = {5(A), Bifi)} = {C(A), C(/x)} = 0, (3.16a) 
{AiX),Bifi)}=-^[Bifi)-BiX)], (3.16b) 

A — fJ, 

{A(A),C(/x)}=-^[C(A)-C(/x)], (3.16c) 

A — H 

{BiX),Ci^^)}=-^[A{^^)-A{X)]. (3.16d) 

A — 

Then {^^(A) + B{X)C{X), A^{fi) + B{fi)C{fi)} = implies that PjJ = 1,...,27V, are 
in involution. The AKNS equation (3.6) is factorized by the x-FDIHS (3.10) and the 
t2-FDIHS (3.11), namely, if (^'i, ^'2, $i, $2) satisfies the FDIHSs (3.10) and (3.11) si- 
multaneously, then {q^r) given by (3.9) solves the AKNS equation (3.6). In general, the 
factorization of the n-th AKNS equations (3.4) will be presented in the end of section 
3.2. 
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3.2 The separation of variables for the AKNS equations . 

In contrast with the B{X) in the Lax matrix M for the constrained KdV flows, the 
B{X) given by (3.12b) has only — 1 zeros, one has to define the canonical variables 
Uk,Vk,k = 1, 2N, in a little different way. The commutator relations (3.16) and the 
generating function of integrals of motion (3.13) enable us to introduce ui,...,un by 
means of -B(A) in the following way: 

«(A) = E|2^^ = e»»|^. (3.17a) 



where 



and vi, Vn by 



,=1 A - A, K{\) ' 



N-l N 

R{\) = J] (A - «fe), K{X) = J] (A - Afe), 



Vk = A{uk), k = l,...,N-l, T;Ar = ^(< > - < *2,*2 >)• (3.17b) 
The equation (3.17a) yields 

= < ^^2, $1 > • (3.17c) 

Then it is easy to verify that 

{uN,B{fj,)} = {vN,A{fj,)} = 0, {vn,un} = 1, (3.18a) 

{UN, A{f,)} = {vn, B{^^)} = 5(/.). (3.18b) 

As we will show later, the commutator relations (3.16) and (3.18) guarantee that 
ui, ...,un, vi, ...jVn satisfy the canonical conditions (1.1). 

We now need to construct two functions B{X), A{X) to define un+i, U2n by means 
of -B(A) and Vjv+i, ■■■,V2n by fAr+fc = A{uN-\-k)- By the exactly same argument as in 
the previous section, we construct A{\) by 

1(A) = \ f ^ f Pn^^ (3.19^) 
2 A — A^ ^ A — A.J 

since the equation (3.19a) enable us to obtain immediately the separated equations (1.2) 
for f jv+fc and UN+k, and -B(A) by 
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Then it is easy to verify that A{X), B{X), A{X), B{pi) satisfy the commutator relations 

{S(A), B{^i)} = {5(A), B{^Ji)} = {1(A), A{^i)} = {A{\),B{^Ji)} = {1(A), A(/.)} = 0, 

(3.20a) 

{A{\),B{^i)} = -^[B{^i) - 5(A)], (3.20b) 

A — )U 

{A(A),5(/.)}= ^[5(/.)-5(A)]. (3.20c) 

The relation (3.20c) doesn't fit the requirement for the canonical conditions (1.1). Ac- 
cording to (3.20) and (3.16) we can replace A{\) by A{\) 

A{X) ^ A{X) - A{X) = -1 - f: (3.21) 

namely, we redefine t^i, vjv by 

vk = A{uk) = A{uk) - Aiuk) = -1 - V ^^HZ^, k=l, N-1, (3.22a) 

vn = - <'i>2,^2> ■ (3.22b) 
We now define un+i, ■■■,U2n by -B(A) as follows: 

^ ^2 _ N-1 

^(^) = E a3T: = '"""ikry ^^^^ = n - ^^+'^)' (3.23a) 
j=l ^ ^ ' k=l 

and vn+1, ■■■,V2n by 

.iv+. = I(^iv+.) = ^E^^^'^^^^^' /c = l,...,Ar-l, (3.23b) 

i^2jv = ^(< > + < *2,$2 >)■ (3.23c) 

The equation (3.23a) leads to 

U2N = In < ^1,^1 > . (3.23d) 

Then a straightforward calculation shows that B{X) = B{X), A{X), B{X), A{X) satisfy 
the commutator relations (2.33) and 

{uN,Bifi)} = {uN,B{i^)} = {uN,A{fi)} = 0, {uN,Aifi)} = -—^i^— , (3.24a) 
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{vN,A{^i)} = {vN,B{fi)} = {vN,A{^i)} = 0, {vN,B{fi)} = Bill), (3.24b) 

{U2N, B{fl)} = {U2N, A{ll)} = {U2N, B{ll)} = 0, {u2N, A{fl)} = > 

(3.24c) 

{v2N,B{ti)} = {v2N,A{fj,)} = {v2N,A{fj,)} = 0, {v2N , B{n)} = B{n) , (3.24d) 

{vn,un}^1, {v2N,U2n} = 1, (3.24e) 

{U2N, Un} = {U2N, Vn} = {V2N, Un} = {V2N, Vn} = 0. (3.24f) 

We have the following proposition. 

Proposition 3. Assume that Xj,(f)ij,i(^ij e R, z = 1,2, j = 1,...,N. Introduce the 
separated variables ui, U2n by the B{X) and B{X): 



N 



S(A)^f:,^=e-»§;. (3.25b) 

N-1 N-1 

R{x) = n - ^fc)' ^(^) = n - ^^+'^)' 

k=l k=l 

and vi, ...,V2N by 

vk = A{uk) = A{uk) - A{uk) = -1 - J] , k=l, N-1, (3.25c) 

J = l 

vjv = - < *2,^2 >, (3.25d) 

vr^+k = Aiu^^k) = \ f: "^'^"^'^ "-^l^f'' , k=l,...,N-l, (3.25e) 



with 



UN+k — \ 



J 

^;2JV = ^(< > + < *2,$2 >)• (3.25f) 

If ui, ...,unj are single zeros of B{X), then vi,...,V2n and ui,...,U2n are canonically 
conjugated, i.e., they satisfy (1.1). 

Proof. By using the exactly same method as in the proof of proposition 1, the com- 
mutator relations (2.33) guarantee that u\, ...,un-i-, v\, ...,vn-i satisfy the canonical 
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conditions (1.1). By using the similar method, for example, it is found from (3.24) that 
for k = 1, N — 1, we have 

= {uN,B{Uk)} = B'{Uk){uN,Uk} + {M7V,5(A^)}|M=Ufe = B' {Uk){uN , Uk} , 

Bifi) , 



{un, Vk} = {un, A{Uk)} = A {Uk){uN, Uk} - 



0, 



which gives rise to {uN,Uk} = {uN,Vk} = and so on. In this way we complete the 
proof. 

It follows from (3.25a) and (3.25b) that 



'fp2j4>ij = e 



or 







— e 



U2N V J 



e"~i?(Aj 



j = l,...,N. 



(3.26) 



The equation (3.9) and (3.17c) results 



q = e 



UN 



(3.27) 



We now present the separated equations. By substituting Uk into (3.13), UN+k into 
(3.19a) and using (3.25c) and (3.25e), one gets the separated equations 



Vk = A{uk) - A{uk) = \/P{uk) - A{uk) 



N 



N 



N+j 



\ ~t^k-^j (Uk-Xj)'^' j^^Uk-Xj 



E^^' ^ = i,...,iv-i, 



N 



VN+k = A{uN+k) = E 



N+j 



^ UN+k — Xj ' 



A;= l,...,Ar- 1. 



(3.28a) 
(3.28b) 



It is easy to see from (3.14) that 



JV 



AT 



i=l 



i=l 



which together with (3.25d) and (3.25f) leads to 



1 ^ 

^^ = -9E(^^ + 2Pjv+0, 



i=l 



N 

V2N = y^^PN+i- 

i=l 



(3.28c) 
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Replacing Vk by the partial derivative of the generating function S of the canonical 
transformation and interpreting the Pi as integration constants, the equations (3.28) 
may be integrated to give the generating function of the canonical transformation 

S{ui,...,U2n)=J2[J iy/P{Xj-A{X))dX + J A{X)dX] 



k=l 

N N 

2 



PN+iU2N 

i=l i=l 
N-1 i,uk N , 

1 -I *^ 1 



^ N N 

- Y,{Pi + 2Pjv+i)«iv + PN+iU2N. (3.29) 



i=l i=l 

The linearizing coordinates are then 

ds 1^ r'^ 1 1 



dS 1 /" 1 1 



(Jjv+i = ttt; = / / , dX — In — — un + W27V, 

^""^^ dPN+i (A-Ai)2/P(A) '^/AT+fc-Ai'J 



i = l,...,iV. (3.30b) 

By using (3.15a), the linear flow induced by the FDIHS (3.10) together with the 
equations (3.30) leads to the Jacobi inversion problem for the FDIHS (3.10) 

N-l uk I N 

J2 / 77 r-j^dX-UN = 'yi + {2Xi-J2Pk)x, i = l,...,7V, (3.31a) 

fe^i ^ {X-Xi)^/P{X) 

N-l 

(A - Xify/P{X)^'^ ' UN+k - \ 

i = l,...,N. (3.31b) 

By using (3.15b), the linear flow induced by the FDIHS (3.11) and the equations 
(3.30) result in the Jacobi inversion problem for the FDIHS (3.11) 

/ dX — um 

tiJ (A-A,)yp(A) 



[ [ — ^y^^-F==«iA - In I \]-Un + U2N = iN+i + 2PN+iX, 



fc=i 
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N N 

= 7i + [2X^i - Y^iXkPk + hPk + PN+k) + ^(JZ Pkf]t2, i = 1, N, (3.32a) 
fc=i fe=i 

\ -^■^+» ^\ ; ~ 1 

/ , aX — In — — UN + U2N 

^1 (A-Ai)2/P(A) UN+k-Xi 



k=l 



N 



^N+i + PN+ii4=Xi - J2 Pk)t2, i = l,...,N. (3.32b) 



fe=i 



Then, since the AKNS equations (3.6) are factorized by the FDIHSs (3.10) and 
(3.11), combining the equations (3.31) with the equations (3.32) gives rise to the Jacobi 
inversion problem for the AKNS equations (3.6) 



^-1 ruk 



^1 (A-A,)v/P(A) 



k=l 

N N ^ N 

= 7i + (2A, - J2 Pk)x + [2X1 - Y^^XkPk + XiPk + P^+fe) + ^($^ Pkf]t2, 

k=l k=l fc=l 



i = l,...,N, (3.33a) 

^ r'^k ID \ 



(A- Ai)2/P(A) UN+k-X 



N 



= ^N+i + 2PN+iX + PN+i{4:Xi-Y,Pk)t2, i=l,...,N. (3.33b) 



fe=i 



If 4>ij,ip2j,<l defined by (3.26) and (3.27) can be solved from (3.33) by using the Jacobi 
inversion technique, then (p2j , ipij can be obtained from the first equation and the last 
equation in (3.10) by an algebraic calculation, respectively. Finally q and r =< ^^i, $2 > 
provides the solution to the AKNS equations (3.6). 

In general, in order to obtain the Jacobi inversion problem for the n-th AKNS equa- 
tions (3.4), we set 

00 

A\X) + S(A)C(A) = J2 FkX-\ (3.34) 

k=0 

where F^-ik = 1,2, are also integrals of motion for both the FDHS (3.10) and the 
tn-binary constrained flow. Comparing (3.34) with (3.13), one gets 

N 

Fo = l, Fk = Yli>^';-'Pj + ik-l)X';-^P^+j], A; =1,2,.... (3.35) 
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The n-th AKNS equations (3.4) are factorized by the a;-FDIHS (3.10) and the t„-FDIHS 
with the Hamiltonian Fn given by [25] 

" 1 ~ ~ 

^n=2 5:(-^)-^ y: ^^.-^w. (3.36) 

m=0 li+...+lm+i=n+l 

where h > 1, Im+i > 1; and am are given by (2.25c). In the same way, we have the 
following proposition. 

Proposition 4. The Jacobi inversion problem for the n-th AKNS equations (3.4) is of 
the form 

N-l u, ^ N 



(A-A,)/P(A) 

+2t„^(--)-a^ \[-^^'-^Fi^...Fi^, z = l,...,iV, (3.37a) 



fe=l " '^t) V ^ k^i 

n 



2' 

m=0 h+...+lm+i=n+l 



(A-Ai)2/P(A) Wiv+fc-Ai 
+4*n^(-2)"«m 5] (Wi-l)A^+^-'Piv+ii^«,...F^^, i = l,...,iV, 

m=0 /i+...+/^+i=n+l 

(3.37b) 

where h > l,...,lm+i > 1, anc? Fi-^,...Fi^, are given by (3.35). 



4. The separation of variables for the Kaup-Newell equations . 
4.1 Binary constrained flows of the Kaup-Newell hierarchy. 

For the Kaup-Newell spectral problem [38] 

</., = [/(«, A) 0, [/(^., A) = '^^(02)' ""^0' ^^'^^ 

its adjoint representation (2.2) and (2.3) yields 

ao = l, a2 = -^qr, bi = -q, Ci = -r, 63 = ^(^^r + g^), C3 = ^(qr^ - r^), 
and in general a2A:+i = ^>2fc = C2fe = 

/C2fc+l^ ^/C2fe-l\ a2fc = ^9"^(5C2fe-l,a: + r62fc-l,x), k = l,2,---, (4.2) 

\02k+lJ \O2k-1J 2 
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^ _ 1 f d — rd~^qd —rd~^rd \ 
2 \ —qd~^qd —d — qd~^rd J ' 

Take 

"iZi / n \2n-2i u \2n-2i-l\ 
i=0 ^ ^"^ * / 



Then the compatibihty condition of equations (4.1) and (4.3) gives rise to the Kaup- 
Newell hierarchy 

[<l\ jfc2n-l\ jSH2n-2 . ^ 

where the Hamiltonian and the Hamiltonian operator J are given by 

j_fQ d\ _ 4a2n+2 - rC2n+l - qhn+l ( C2n+l\ _ ^-^2n 

For 77. = 2 we have 

and the coupled derivative nonlinear Schrodinger (CDNS) equations obtained from the 
equation (4.4) for n = 2 read 

(lt2 = \lxx + ]^{.Q^r)j;, rt^ = ~rxx + ^{r'^q)x. (4.6) 

The adjoint Kaup-Newell spectral problem is the equation (2.11) with U given by 
(4.1). We have [25] 



^ f^Trlf'f'^'f'^ 01^2 \ ^^(^,A) _ /A^102\ . . 

The binary x-constrained flows of the Kaup-Newell hierarchy (4.4) are defined by 

$1,0; = -A^^i +gA$2, $2,0; = rA$i + A2$2, (4.8a) 
= A^*! - rA*2, *2,a; = -gA*i - (4.8b) 



N 



SHko SXj fc2ko + l\ 1 / < A^'i, $2 > 
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For A;o = 1, we have 

By substituting (4.9) into (4.8a) and (4.8b), the first binary a;-constrained flow becomes 
aFDHS 

dFi ^ dFi ^ dFi ^ dFi 

1^ 5^1 d^2 

with the Hamiltonian 

Fi =< A^*2,^2 > - < A^*i,$i > < A*2,^i >< A*i,^2 > • 

Under the constraint (4.9) and the FDHS (4.10), the binary t2-constrained flow ob- 
tained from (4.3) with V^'^'^ given by (4.5) and its adjoint equation for N distinct reral 
numbers A,- can also be written FDHS 

dF2 ^ dF2 ^ dF2 ^ dF2 , 

with the Hamiltonian 

I 

F2 = - < A^*2, $2 > + < A^*i, $1 > +7: < A*2, $1 >< A^*i, $2 > 
1 1 

+ - < A^*2, ^1 >< A*i, $2 > < A*2, ^1 >^< A*i, $2 >^ 

I 

+ -(< A^*2, $2 > - < A^*i, $1 >) < A*2, ^"1 >< A*i, $2 > • 

8 

The Lax representation for the FDHSs (4.10) and (4.11) are presented by (2.18) with 
the entries of the Lax matrix M given by 

1 ^ 1 V" ^^(V'lj^ij - V'2i</>2i) 
BW = iAf:M^, C(A) = iAf:^|^. (4.12b) 

j=l j j=l 3 

A straightforward calculation yields 

N p A'^P^ 
^^(A) + 5(A)C(A) ^ P(A) = 1 + E[a^ + {X^-X]r ^' ^^'^^^ 
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where Pj,j = 1, ...,2N, are 2N independent integrals of motion for the FDIHSs (4.10) 
and (4.11) 



+ J TJ3^[AjA|(V'lj01j-V'2j02jO(V'lfc01fc-V'2fe02fc) + 2AjA^^ 

j = l,...,N (4.14a) 

Pn+j = + V'2,</'2j), j = 1, N. (4.14b) 
It is easy to varify that 

TV N ^ N 

F, = -2j2Pv F2 = 2^(A,2p,. + A|p2^,.) - -(5]P,)^ (4.15a) 

AT 

< ^'2,^2 > + < >= 4 J^Ptv+j. (4.15b) 
By inserting A = 0, (4.13) leads to 



1 + ^(< *2,$2 > - < >) = /P(0) = 



TV 



With respect to the standard Poisson bracket it is found that 

{A{\),A{^Ji)} = {BiX), Bifi)} = {C(A), Ci^^)} = 0, (4.17a) 

{A{X), B{^)] ^ [^B{^) - AP(A)], (4.17b) 

{A(A), C(^)} = 2(A2^1^2) [^^(-^) - ^^(a*)]' (4-17c) 

{P(A),C(/.)} = ^^5^[^(^) -^(A)]. (4.17d) 

Then {^^(A) + 5(A)C(A), ^^(/x) + B{tJ,)C{fj,)} = implies that PjJ = 1, ...,2iV, are 
in involution. The CDNS equations (4.6) are factorized by the x-FDIHS (4.10) and 
the t2-FDIHS (4.11), namely, if (^'i, ^'2, $2) satisfies the FDIHSs (4.10) and (4.11) 
simultaneously, then (g, r) given by (4.9) solves the CDNS equations (4.6). The factor- 
ization of the n-th Kaup-Newell ewuations (4.4) will be presented in the end of section 
4.2. 
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4.2 The separation of variables for the Kaup-Newell equations. 



Since the commutator relations (4.17) are quite different from (2.23) and (3.16), we 
have to modify a Httle bit of the method presented in sections 2 and 3. Let us denote 
X — X^, Xj — Xj. The entries of the Lax matrix M given by (4.12) can be rewritten as 

A(X) = 1 + |(< ^2, $2 > - < ^-1,$! >) + ^AAi(A), b(X) = ^Vxb(X), (4.18a) 



where 



N 



1 ^ VlyC'l.v - V', 



N 



A- A, 



j=i A - Xj 



It is easy to see that 



{A^iX),A^i^l)} = {B{X),Bi|l)} = 0, 

{A,cx),m} = r^[m-B{x)]- 



X — jJL 



(4.18b) 

(4.19a) 
(4.19b) 



It foUows from (4.16) and (4.18a) that 
A{X) 



\ 



(4.19c) 



i=i 



The commutator relations (4.19) and the generating function of integrals of motion 
(4.13) enable us to introduce ui, ...,un in the following way: 



N 



^j'4'2j4>ij 



,=1 A - Aj 



K{X) 



(4.20) 



where 



N-l 



N 



fe=i 



fe=i 



and Vat by Ai(A). 

By the exactly same argument as in sections 2 and 3, we construct A{X) and B{X) 

by 



N 



N 



A-A, fctA-A,' 



B(A) = E 



,=1 A - Xj 



(4.21) 
(4.22) 
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and, for the same reason, we have to replace Ai{X) by A{X) 

A(A).A,(A)-I(A) = -f:|^. (4.23) 
Then we have the following proposition. 

Proposition 5. Assume that \j^4'iji'4'ij ^ F^? * = l)2,j = 1,...,N. Introduce the 
separated variables ui, U2n by the B(X) and B{X): 



j=i A - K (A) 

BiX) = y = (4.24b) 
,ttA-A, K{X) 

with 

N-l N-l 

r(X) =l[(X-Uk), R(^) = n - 

k=l k=l 

and vi, ...,V2N by 

Vk = A{uk) = A,{uk) - A{uk) = -J2 k = l, N-l, (4.24c) 

vtv = - < *2,$2 >, (4.24d) 



= !(..+,) = ^f:^^^M^i^±%^, k=l,...,N-l, (4.24e) 

j=l UN+k — 

V2N = ^(< *!,$!> + < *2, $2 >)■ (4.24f) 

If ui, ...,un, are single zeros of B{X), then vi,...,V2n and ui,...,U2n are canonically 
conjugated, i.e., they satisfy (1.1). 

Proof. It follows from (4.24a) and (4.24b) that 

UN = ln< A*2, *i >, (4.25) 

U2N = < *i, $1 > . (4.26) 
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By a straightforward calculatfon, it is found that -B( A) , ^(A) , -B( A) , A{X) satisfy the com- 
mutator relations (2.33) with A, replaced by A, Jl, as well as the following commutator 
relations 

{uN, Bill)} = {UN, Bill)} = {uN,A{fi)} = 0, {uN,Aiii)}= ^^^^ 



< A^'2,$l >' 

(4.27a) 

{vN^Aii,)} = {vN,Biij,)} = {vN^Aii,)} = 0, {vN,Bii,)} = Bill), (4.27b) 

{u2N,Bili)} = {u2N,Ain)} = {U2N, -B(/i)} = 0, {U2N, Mli)} - ^^^^ 



(4.27c) 

{v2N.Bii,)} = {v2N,Aiii)} = {v2N,AiiJ,)} = 0, {V2N . B ill)} = B ill) , (4.27d) 

{vn,un} = 1, {v2n,U2n} = 1, (4.27e) 

{U2N, Un} = {U2N, Vn} = {v2N, Un} = {v2N, Vn} = 0. (4.27f) 

Then in the exactly same way as for the propositin 1 and 3, we can complete the 
proof. 

It follows from (4.24a) and (4.24b) that 



i?(A2) „ , i?(A2) 

^ J,2 _ U2N ^ J ' 



XjiP2jcf>ij = e"-— = e'^^'T^, j = 1, -.N, (4.28) 



or 



'e"2^i?(A2) _ e"^i?(A2) 

^'(^^) ' Xj^e^^^RiX])K'iX]) 



^^J = \ ^/M2^ > ; _ ' j = l,...,N. (4.29) 



The equations (4.9) and (4.25) result 

g = -le"^. (4.30) 

We now present the separated equations. By substituting uj. into (4.13), UN+k into 
(4.21) and using (4.19c), (4.24c) and (4.24e), one gets the separated equations 

Vk = A^iuk) - Aiuk) = —[\IPink) - \fPm - -^I^' ^ = 1' -'^ - 1' 



Uk ^ ^,Uk-X' 



(4.31a) 



N 



VN+k = AiuN+k) = 2 V k=l, TV - 1, (4.31b) 
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where -P(O) are given by (4.16) and 

It follows from (4.15b), (4.16), (4.24d) and (4.24f) that 

TV AT 

^;jv = 2-2v/p(0)-2j]Pjv+i, V2JV = 2 J] Pjv+i- (4.31c) 

Replacing by the partial derivative 7^ of the generating function 5" of the canonical 
transformation and interpreting the Pi as integration constants, the equations (4.31) 
may be integrated to give the generating function of the canonical transformation 

S{ui,...,U2n)^J2U (-V^(A)--/P(0)-^(A))o?A+ / A{X)d\] 

A;=l A A J 



N N 



+ (2-2yP(0)-2 5^ PN+i)uN + 2 ^ PN+iU2N 
i=l i=l 

N-1 .... JV , 



^^PiX)dX-2y/P(^ln\uk\-2j2PN+iln 
A 



^2 

k=l " 1=1 --.v-rn, • ^ 



UN+k — \ 

+ (2 - 2/P(0) - 2 ^ Pjv+i)«iv + 2 Piv+i«2Ar. (4.32) 



AT A? 



The linearizing coordinates are then 

dS ^\ 1 ~ 1 , , „ 1 



ti^ A(A-Af)V^ 



UN, 



i=l,...,N, (4.33a) 

^ p 2AfP,v+, ^~ 2P^ , u. A? . 



Pn+i 



-2{-^^ + l)uN + 2u2N, i = l,...,N. (4.33b) 
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By using (4.15a), the linear flow induced by (4.10) together with the equation (4.33) 
leads to the Jacobi inversion problem for the FDIHS (4.10) 



[ / dX H /nl^fcl] -\ UN = li — 2x, 



i=l,...,7V, (4.34a) 



f / — ^! d\ , '^^ ln\uk\\ — 2ln , 



-2{^^ + l)uN + 2u2N = lN+i, i = l,...,N. (4.34b) 

By using (4.15a), the linear flow induced by (4.11) and the equation (4.34) yield the 
Jacobi inversion problem for the FDIHS (4.11) 

1 ~ 1 . , „ 1 



fc=W A(A-Af)^P(A) \VPW KVPW k=i 

i=l,...,N, (4.35a) 

fcli A(A-Af)2yp(A) UN+k-X, 

-2{^^ + l)uN + 2u2N=7N+i+^XjPN+it2, i = l,...,N. (4.35b) 

Finally, since the CDNS equations (4.6) are factorized by the FDIHS (4.10) and 

(4.11), combining the equation (4.34) with the equation (4.35) gives rise to the Jacobi 
inversion problem for the CDNS equations (4.6) 



^-'^ ruk 1—1 1 

/ [/ , dX -\ — - — , ^nlrtfcl] H — - — , UN 

N 

= 7i-2x+(2A^- J]Pfe)t2, i=l,...,N, (4.36a) 



fc=i 

N-l 



E [ r dX - ^^ln\u,\] - 2ln \ ~ 

tiJ XiX-Xir^) ^iv+.-Af 
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-2{^^ + l)uN + 2u2N^lN+i+^XjPN+it2, i = l,...,N. (4.36b) 

If (j)ij,ip2j,q defined by (4.29) and (4.30) can be solved from (4.36) by using the Jacobi 
inversion technique, then 02j,V'ij can be obtained from the first equation and the last 
equation in (4.10), respectively. Finally q and r = — < A^^i, $2 > provides the solution 
to the CDNS equations (4.6). 

In general, the above procedure can be applied to the whole Kaup-Newell hierarchy 
(4.4). Set 

00 

A\\) + B{\)C{\) = J2Fk)^~^\ (4.37a) 

fc=0 

where Fk, k = 1,2, are also integrals of motion for both the x-FDHSs (4.10) and the 
tn-binary constrained flows (2.16). Comparing (4.37a) with (4.13), one gets 

N 

Fo = l, Fk = J2[Xf-^Pj + {k-l)XfP^^+^], A; =1,2,.... (4.37b) 

i=i 

By employing the method in [34,35], the tn-FDIHS obtained from the tn-constrained 
flow is of the form 

with the Hamiltonian 

n — l ^ 

^n = 2 5](-i)--^ Y: (4-38b) 

m=0 ii + ...+Z^+i=n 

where li > 1, ...,lm+i > 1, and are given by (2.25c). Since the n-th Kaup-Newell 
equations (4.4) is factorized by the x-FDIHS (4.10) and the in-FDIHS (4.38). We have 
the following proposition. 

Proposition 6. The Jacobi inversion problem for the n-th Kaup-Newell equations (4-4) 
is given by 

1 ,r 1 , , . 1 



[ / dX H ^nl^fcl] H um = 7i — 2x 

n—l ^ 

+2*n^(-2r«m Yl Xf-+'-^Fi,...Fi^, i=l,...,N, (4.39a) 

m=0 li + ...+lm+i=n 
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-2( ^^^^ + l)«iv + 2u2iv = 7iv+i 

n—l ^ 

+^tr,J2{--ram Yl {lm+i-l)\f"'^'PN+rFi,...Fi^, ^=l,...,iV, (4.39b) 
m=0 h + .-.+^m+i^n 

where /i > 1, > 1? ond F^j, ■■■Fi^, are given by (4- 37b). 

For example, the third equations in the Kaup-Newell hierarchy with n = 3 are of the 
form 

13 13 

qts = -^Qxxx - -^{q^r'^ + '2qrqx)x, = --r^xx - -^{r^Q^ - '2qrrx)x- (4.40) 

The Kaup-Newell equations (4.40) can be factorized by the x-FDIHS (4.10) and ^3- 
FDIHS with the Hamiltonian F3 defined by 

N N N N 

Fa = $:(2A|P, + 4A,^P^^,.) - [^^(A^P, + A|P^+,.)] E + ^(E P^)'- (4.41) 
The Jacobi inversion problem for the equations (4.40) is given by 

N-l "Uk T _ I 

dX + ——==ln\uk\] + ——==un 



/I 
A(A-A^)V^^'" ' Afym^'^'^'^'^ ■ Afv^ 



AT N 

7, - 2a; + [2Xt - Y.^X]P^ + + + ^ (E ^i)']^3, ^ = 1, 

E [ r ~ ~ ^^^^"^^^ ^^ - '^^HukW - 2ln ' "'^"^^ 



fe=i"^ A(A-Af)2W^ UN+k-Xi 



P ^ 
-2(— + 1)UN + 2«2Ar = 7iv+i + [SAfPiv+i - 2A^PAr+j E ^j]*3, ^ = 1, N. 
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5. Concluding remarks. 

The method in [2,3,4,8] is not vahd for the separation of variables for binary con- 
strained flows of sohton hierarchies. This paper proposed a new method to solve this 
problem. For a certain kind of integrable models, a general approach to the separa- 
tion of variables was proposed in [2,3,4] by taking the poles of the properly normalized 
Baker- Akhiezer function and the corresponding eigenvalues of the Lax operator as sepa- 
rated variables. As pointed out in [4] , there is no guarantee that the separated variables 
so constructed satisfy the canonical conditions (1.1). The method proposed in this pa- 
per is to start directly from the canonical conditions (1.1) and the requirement for the 
separated equations (1.2). We introduced 2N pairs of separated variables by means 
of four functions B{X), A{\),B{\) and A(A), which are constructed in such way that 
they satisfy certain commutator relations required by the canonical conditions (1.1) and 
A{\) and ^(A) are linked to the generating functions of the integrals of motion for the 
models. This method ensures that the separated variables are canonically conjugated. 
We produced two sets of separated equations directly from two generating functions of 
intergrals of motion. It seems that the separated equations are intimately connected 
with the generating functions of intergrals of motion. 

The finite gap solutions or finite-dimensional quasiperiodic solutions for the KdV 
equation was studied in [36] by means of the stationary equation of the KdV hierarchy 
called the Lax-Novikov equation. By the standard Jacobi inversion technique [19], the 
finite-dimensional quasiperiodic solution can be given in an explicit form in terms of the 
Riemann theta functions associated with the invariant spectral curve. The Jacobi in- 
version problem (2.45) for the KdV equation is somewhat different from that in [36] due 
to some additional terms. The Jacobi inversion problems for some binary constrained 
flows require the development of the standard Jacobi inversion technique in order to 
solve them explicitly. 
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